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Based on the well-known relation between Fokker-Planck equations and Schrodinger equations of 
quantum mechanics (QM) , we propose new deformed Fokker-Planck (FP) equations associated with 
the Schrodinger equations of "discrete" QM. The latter is a natural discretization of QM and its 
Schrodinger equations are difference instead of differential equations. Exactly solvable FP equations 
are obtained corresponding to exactly solvable "discrete" QM, whose eigenfunctions include various 
deformations of the classical orthogonal polynomials. 

PACS numbers: 05.40.-a, 02.30.1k, 03.65.-w , 02.30.Gp 



-a 
c 

o 
o 



X 



1. Since it was first introduced by Fokker and Planck 
to describe the Brownian motion of particles, the Fokker- 
Planck (FP) equation has become one of the basic tools 
used to deal with fluctuations in various kinds of systems 
[lj. Recently, the phenomena of anomalous diffusion in 
fractal and disordered media have prompted new devel- 
opments in the FP theory. Most attempts have been 
along developing fractional FP equations, where the ordi- 
nary spatial derivatives are replaced by fractional deriva- 
tives Q. Others have tried to generalize the linear FP 
equation to nonlinear ones Q . Diffusion equation based 
on ^-derivatives has also been considered [4J. 

In this paper, we shall generalize the FP equation 
in a different direction. We shall derive new types of 
deformed FP equation which are associated with the 
Schrodinger equations of the "discrete" quantum me- 
chanics considered in @, @ . The eigenfunctions of some 
exactly solvable "discrete" quantum mechanics include 
various deformations of the classical orthogonal poly- 
nomials (the Hcrmitc, Laguerre and Jacobi polynomi- 
als), namely, those belonging to the family of Askey- 
scheme of hypergeometric orthogonal polynomials, and 
the Askey- Wilson polynomial in g-analysis @, H, @ ■ The 
discrete Schrodinger equation is a discrete deformation 
of the usual Schrodinger equation in much the same way 
that the Ruijsenaars-Schneider-van Diejen type systems 
[l(| are the "discrete" counterparts of the Calogero and 
Sutherland systems [TTj ] . the celebrated exactly solvable 
multi-particle dynamics. In fact, these deformed orthog- 
onal polynomials arise in the problem of describing the 
equilibrium positions of Ruijsenaars-Schneider type sys- 
tems @, corresponding to the facts that the Hermite, 
Laguerre and Jacobi polynomials describe the equilib- 
rium positions of the Calogero and Sutherland systems 
In this sense, our new deformed FP equations can 
be considered as the "discrete" deformation of the usual 
FP equation. 



2. In one dimension, the FP equation of the probability 



density P(x,t) is [l[ 



dP(x,t) 



Of 

L FP = -£-DW(x) 



L FP P(x,t), 



dx 
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The functions D^(x) and D^ 2 \x) in the FP operator 
L F p are, respectively, the drift and the diffusion coef- 
ficient (we consider only time-independent case). The 
drift coefficient represents the external force acting on 
the particle, while the diffusion coefficient accounts for 
the effect of fluctuation. The drift coefficient is usually 
expressed in terms of a drift potential <$>(x) according to 
D {1 \x) = -$'(x), where the prime denotes derivative 
with respect to x. 

It is well known that the FP equation is closely re- 
lated to the Schrodinger equation. The correspondence 
between the two equations is usually given by transform- 
ing the FP equation into the corresponding Schrodinger 
equation [l| . But the argument could well be reversed. 

Consider a quantum mechanical Hamiltonian of one 
degree of freedom (we adopt the unit system in which h 
and the mass m of the particle are such that h = 2m = 1) 



H=-^ + V [x) . 



(2) 



Let V{x) be such that the ground state energy is zero, i.e. 
Hcj)^ — 0. The Hamiltonian H is completely determined 
by its ground state wave function 4>q(x). By the well- 
known theorem of quantum mechanics, 4>o{x) has no node 
and can be chosen real. Thus it can be parametrized by 
a real prepotential W(x): 



4> (x) 



V{x) = W'{x) 2 -W"(x). 



(3) 
(4) 



The Hamiltonian is factorizable and positive semi- 
definite: H = A^A. Here A^ = -d x + W and A = 
d x + W. The constant part of W(x) is so chosen as to 
normalize 4>o(x) properly, J (j)o(x) 2 dx — 1. 

Now we define an operator from the H and </>o by the 
similarity transformation 



L 



FP 



-<f>oH(j> 1 , 



(5) 



2 



which guarantees the non-positivity of the eigenvalues of 
Lpp. With (f>o given in ([3]), one obtains 



d , d 2 
L FP = —2W' + ^. 
ox Ox z 



(0) 



From ((TJ), it is seen that Lfp is just the corresponding FP 
operator with = 2W as the drift potential and with 
a constant diffusion coefficient (here equals one owing to 
the choice of unit in H). Hence both H and Lpp are 
determined by 4>o- The eigenfunction P n (x) of Lpp cor- 
responding to eigenvalue — A„ is related to the (real and 
normalized) eigenfunction 0„ of H corresponding to A„ 
by P n {x) = 4>o(x)4> n (x). The stationary distribution is 
Pa = 4>q = exp(— 2W), which is obviously non-negative, 
and is the zero mode of Lpp: LppPo = 0. Any pos- 
itive definite initial probability density P(x, 0) can be 
expanded as P(x 1 0) = 4>o(x)^2 n c n 4>n(x), with constant 
coefficients c n (n = 0, 1, . . .) 



0„(x) {^ 1 {x)P{x,Q))dx. 



Then at any later time t, the solution of the FP equation 

is P(x,t) = 4>f){x)Y, n Cn4>n{x) 6Xp(-A n t). 

3. We now derive a new class of deformed FP equa- 
tion corresponding to the discrete Schrodinger equations 
discussed in [f| in accordance with the prescription in 
the last section. Eigenfunctions in this type of discrete 
Schrodinger equations are related to the family of Aske^ 
scheme of hypergeometric orthogonal polynomials @, \ 

The Hamiltonian has the general form 



H 



jV&j e~ w * ^/VHx) + VV^ix) e w * ^V(x) 
-(V(x)+V*(x)). 



(8) 



Here the momentum operator p = —id x (with K = 1) 
appears as exponentiated instead of powers in ordinary 
quantum mechanics. Thus they cause a finite shift of the 
wave function in the imaginary direction: e^ 19 * (f>(x) — 
4>(x±i). Throughout this paper the following convention 
of a complex conjugate function will be adopted: for an 
arbitrary function f(x) — ~Y^ n a n x n , a n <E C, we define 
f*( x ) — En "n 1 "- Here a* is the complex conjugation 
of a n . Note that f*(x) is not the complex conjugation of 
f(x), {f(x))* = f*(x*). This is particularly important 
when a function is shifted in the imaginary direction. 
The Hamiltonian 15]) is factorised, i.e. H = A^A with 
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y/V*{x) 



e 2 



VV*{x), 



(9) 
(10) 



Here f denotes the ordinary hermitian conjugation with 
respect to the ordinary L 2 inner product: (f\g) — 
J^^ifix))* g(x)dx. Obviously the Hamiltonian © is 
hermitian (self-conjugate) and positive semi-definite. 
Let the ground state 0o be annihilated by A: 

A(j> {x) = 0, (=► HM*) = 0, A = 0). (11) 



Explicitly the above equation reads 



.JV* {x - i/2) O (x - i/2) = ^/V(x + i/2) O (x + i/2) , 

(12) 

or equivalently 



^V*{x-i) 4>o{x~i) = y/V{x)Mx), (13) 
y/V(x + i) 4>o{x + i) = ^V*{x) o (z). (14) 

Eq. (fT2")) relates the potential V(x) and the ground state 
0o, and hence is the discrete analogue of (J3j) albeit in an 
implicit way. 

Now, we form the associated FP operator from ([SJ and 
0o in (fTTj) according to the similarity transformation |(SJ). 
The sum V(x) + V*(x) remains intact. The first term in 
H transforms as 

0o (x) y/V{x) e~ id ° ^/W{x) O - 1 (a;) 



<h{x-i)y/V*{x-i) 



V*(x). (15) 



By virtue of (flU)) . the term in the bracket of the r.h.s. of 
(I15p is simply unity. Hence 



0o (x) ^/V{x) e~ w * VV*(x) 0o 1 (x) 



-id. 



V*(x). (16) 



Similarly, using (|14|l , the second term in H transforms as 
Mx)^VHxje w *^V(x) o ~ 1 (x) = e id *V(x). (17) 
Putting everything together, the result is 

Lpp = -e l9 *V(x) - e- i9 *V*{x) + V(x) + V*(x). (18) 

This is the general form of FP operator corresponding to 
the discrete Hamiltonian H in |8j). One has Lpp§\ = 
as a consequence of i?0o = 0. Thus 0q is the stationary 
solution of the respective FP equation. 

We now discuss the limiting form of the FP equation 
when the momentum p = —id x is small. We expand the 
operators e ±ld:c in (fT5|) up to the 2nd order in p. The FP 
operator becomes 

Lpp = -r^Im V(x) + g^ZRe V(x). (19) 

Hence, in this limit, the deformed FP equation (fT8"|) 
does reduce to the usual FP equation {T]), with D^> — 
-21m V(x) and = 2Re V{x). 

Let us illustrate this connection with a simple example 
discussed in [5j], namely, the discrete Schrodinger equa- 
tion with the Meixner-Pollaczek polynomials as eigen- 
functions. These polynomials are deformation of the 
Hermite polynomials. The potential in this system is 
V(x) = a + ix with a real and positive. From the above 
discussion, the small momentum limit of the FP equation 
is just the ordinary FP equation with = —2x and 
Z)( 2 ) = 2a. This latter system is none other than the FP 
equation of the celebrated Ornstein-Uhlenbeck process, 
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which can be exactly solved by means of expansion in 
terms of the Hermite polynomials [l[ . 

4. In Q four shape-invariant potentials [l3| in the dis- 
crete quantum mechanics were presented, including the 
Meixner-Pollaczek case discussed above. Shape invari- 
ance ensures that the corresponding Schrodinger equa- 
tions are exactly solvable. Hence, the corresponding 
deformed FP equations are also exactly solvable. For 
completeness, we list below the potentials V(x), the cor- 
responding ground state wave functions (f>o(x), and the 
eigenvalues A„ (n — 0, 1, 2, . . .). Needless to say all the 
eigenf unctions are square integrable. All parameters a, 
b, c and d are assumed to be real and positive. We name 
the case by the type of polynomials to which the poly- 
nomial part of 4> n belongs (see [|| for the details of these 
polynomials). 

(i) the Meixner-Pollaczek case, 

(one-parameter deformation of the Hermite polynomial) 

V(x) = a + ix, 
(j)o(x) oc \T(a + ix)\, 
A„ = 2n; 

(ii) the continuous Hahn case, 

(two-parameter deformation of the Hermite polynomial) 

V(x) = (a + ix)(b + ix), 
(j)o(x) oc \T(a + ix)T(b + ix)\, 
X n = n(n + 2a + 2b - 1); 

(iii) the continuous dual Hahn case, 
(two-parameter deformation of the Laguerre polynomial) 



V(x) 
<j>o(x) oc 



(a + ix){b + ix)(c + ix) 
~ 2ix{2ix + 1) ' 

r(a + ix)T{b + ix)T(c + ix) 



r(2ix) 



A„ = n; 



(iv) the Wilson case, 

(three-parameter deformation of the Laguerre polyno- 
mial) 



V(x) = 



(a + ix)(b + ix)(c + ix)(d + ix) 



<po(x) oc 



2ix(2ix + l) 
T(a + ix)T(b + ix)T(c + ix)T(d + ix) 



T(2ix) 

X n = n(n + a + b + c + d— 1). 



For compactness, we have written the ground state wave 
function in terms of the absolute value symbol. For in- 
stance, 



\T(a + ix)\ = yT(a + ix)T(a — ix). 



5. The potentials listed in the last section are exactly 
solvable as they are shape-invariant. However, the con- 
struction of the deformed FP operator can be carried over 
to the more general potential 



V(x) = 



ILK 



Y\ k {d k +ix)' 



(21) 



where dj, d k are real and positive, and the degree of 
the numerator is at least one greater than that of the 
denominator for square integrability. The corresponding 
ground state wave function is 



(f> (x) cx 



/ Ylj T(cij + ix)T(aj — ix) 
Yl k T{d k +ix)T(d k -ixY 



(22) 



The potential V(x) and the state (f>o(x) satisfy the rela- 
tion (fP2"|) . As discussed previously, this ensures that the 
Hamiltonian H with V(x) in (|2ip is mapped into Lfp in 
(JT5J) with the similarity transformation ^ , and that the 
probability density <pQ is indeed the zero mode of Lfp- 

6. We turn now to a different type of shape-invariant 
"discrete" quantum mechanical single particle systems 
with a g-shift type kinetic term discussed in Q. Eigen- 
functions in this class of system are related to the Askey- 
Wilson polynomials (sj. We now derive the associated 
FP operator for such systems. 

Following Q, we use variables 8, x and z, which are 
related as 



< < it, x = cos6, z = e l 



(23) 



The dynamical variable is 9 and the inner product is 
= f*ddf(&)*g(6). We denote D = z£. Then q D 
is a g-shift operator, q D f(z) — f(qz). We note here that 



dO = 



dx 



-iVT 



- l Te = z i = D ' f^* = f*^ 1 )- ( 24 ) 

For a function V(z), which is a function of a real con- 
stant q (0 < q < 1) and a set of real parameters, we 
define the following Hamiltonian H, 



(20) 



H = ^V{z)q D ^V*(z-^ + ^/V*{z^) q- D ^V{z) 

-(V(z) + V*(z- 1 )). (25) 

The eigenvalue equation reads H<j) n — \ n 4> n with eigen- 
functions 4> n (z) and eigenvalues A„ (n = 0,1,...) (we 
assume non-degeneracy Ao < Ai < •■•). The kinetic 
term causes a q-shift in the variable z. This Hamiltonian 
is factorized, i.e. H = A^A, with 

A = q^^V^z-^-q-^^Vir), (26) 
At = ^V\z)q% - ^V^z-^q-^. (27) 
The ground state <po is the function annihilated by A: 
A(j) = (=► Hfo = 0, A = 0). (28) 
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Explicitly this equation reads 

yfviq-tz-^Ml**) = \fv(q-tz)Mq-*z)- (29) 
Note that Eq. ([29]) implies 

^V*{q-^Z-^) O (qz) = y/V{z) (30) 
y/V*(z-i) <h{z) = y/Viq-iz) <M<rH (31) 

As with fl2]l. ([29]) is the analogue of ©. The other 
eigenfunctions can be obtained in the form 



4> n {z) oc p n (z)0o(z), 



(32) 



where p n (z) is a Laurent polynomial in z 

Using iJ, 0o, ([30)) and ([3~T]) . the similarity transforma- 
tion © produces the FP operator 



L 



FP 



-D 



V(z) — q D V*(z~ 1 ) 



+V{z) + V*{z' 1 ). 



(33) 



This is the general discrete cy-deformed FP operator cor- 
responding to the Hamiltonian (|25|) . Again, Lpp annihi- 



lates 0o- 



As an example, let us take V(x) to be the one discussed 



= {l-az){l-bz){l-cz){l-dz) 
(1 — z z )(l — qz z ) 

For simplicity we assume — 1 < a, b, c, d < 1. Note that 
V*(z) = V(z). The ground state is given by [|[ 



0o (z) oc 



(* 2 ;</)c 



(az, 6z, cz, dz; g) c 



(35) 



(z 2 ,z 2 ;q) 



(az,az 1 ,bz,bz 1 ,cz,cz 1 ,dz,dz 1 ',q) c 



where (ai,--- ,fl m ;9) M = n, = i Yl n =o( 1 ~ a j^ n )- Ex cited 
states have the form ([32]) <f) n (z) oc p n (z)4>o(z), where 
Pn(z) is proportional to the Askey- Wilson polynomial 
[g, [9J , which is a three-parameter deformation of the Ja- 
cobi polynomial. The eigen-energies are 



q- n (l - q n )(l - abedq 71 ' 1 ) 



(36) 



Being shape- invariant, the Schrodingcr equation of this 
system is exactly solvable, and so is the corresponding 
FP equation. 
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